Effects of Strong Color Fields on Baryon Dynamics 
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We calculate the antibaryon-to-baryon ratios, p/p, A/A, H/H, 
and U/D. for Au+Au collisions at RHIC (Vsjvjv = 200 GeV). 
The effects of strong color fields associated with an en- 
hanced strangeness and diquark production probability and 
with an effective decrease of formation times are investi- 
gated. Antibaryon-to-baryon ratios increase with the color 
field strength. The ratios also increase with the strangeness 
content The netbaryon number at midrapidity consid- 
erably increases with the color field strength while the net- 
proton number remains roughly the same. This shows that 
the enhanced baryon transport involves a conversion into the 
hyperon sector {hyperonization) which can be observed in the 
(A — A)/(p — p) ratio. 

Ultrarelativistic heavy ion collisions provide a unique 
tool to study elementary matter at energy densities so 
high that a phase transition from partonic deconfined 
matter to hadronic matter is predicted by QCD lattice 
calculations. The properties of such a highly excited 
state depend strongly on the initial conditions such as 
the specific entropy density or the thermalization time. 
A quantum number of utmost interest is the baryon num- 
ber which, prior to the collision, is non-vanishing only 
at beam and target rapidities which for energies at the 
Relativistic Heavy Ion Collider (RHIC) (center-of-mass 
energy per nucleon-nucleon collision -v/sAfw ^ 

200 GeV, 

Au-|-Au) are separated by 

Ay ~ 11 

units. Accordingly, one often assumes complete trans- 
parency, i.e., vanishing netbaryon number and B/B- 
ratios identical with one at midrapidity. Recent data by 
the STAR [1], PHOBOS [2], BRAHMS [3], and PHENIX 
[4] collaborations, however, show that the midrapidity re- 
gion 2/^0 obtains a considerable finite netproton {p -— p) 
number in the final state. Moreover, the antibaryon-to- 
baryon ratios are not yet equal to one. The p/p-ratio is 
about 0.7 for the most central collisions. Furthermore, 
rapidity spectra seem to exhibit rather a plateau than 
a minimum at midrapidity and large elliptic flow values 
indicate the need for a large amount of momentum equili- 
bration at midrapidity providing additional evidence for 
a non-transparent scenario. 

This strong baryon number transport over more than 
5 units in rapidity is a novel phenomenon and its under- 
lying mechanism is of utmost interest [5-9] . It also has a 



strong impact on other questions like the (parton) equili- 
bration [10], the equilibration in terms of hadronic yields 
[11] and the build up of collective flow [12]. 

In this Letter we study the effects of strong color fields 
on the dynamics of baryons in relativistic heavy ion col- 
lisions. We will show that strong color fields lead to 
an enhanced baryon-antibaryon pair production and to 
a strong modification of the baryon transport dynam- 
ics. This will affect the antibaryon-to-baryon ratios as 
well as the netproton and netbaryon rapidity distribu- 
tions. Moreover, the netbaryon transport to midrapidity 
includes a hyperonization, i.e., the conversion of ordinary 
nonstrange baryons to baryons carrying strangeness. 

The idea of strong color electric fields, i.e., of an ef- 
fectively increased string tension in a densely populated 
colored environment of highly excited matter has been 
suggested already earlier. The particle formation process 
in hadronic collisions can be viewed as quantum tunnel- 
ing of quark-antiquark and gluon pairs in the presence 
of a background color electric field. It is formed between 
two receding hadrons which are color charged by the ex- 
change of soft gluons while colliding. In nucleus- nucleus 
collisions the color charges may be considerably greater 
than in nucleon-nucleon collisions due to the almost si- 
multaneous interaction of several participating nucleons 
[13]. This leads to the formation of strong color electric 
fields. With increasing energy of the target and pro- 
jectile the number and density of strings grows, so that 
they start overlapping, forming clusters, which act as new 
effective sources for particle production [14,15]. It has 
been predicted that the multiplicities of, for example, 
strange baryons or antibaryons should be strongly en- 
hanced [15-18] once the color field strength grows. The 
abundances of (multiply) strange (anti)baryons in central 
Pb-|-Pb collisions at Cern-SPS [19], for example, can only 
be explained within the framework of microscopic model 
calculations [18] if the elementary production probabil- 
ity of ss pairs, which is governed in the string models 
[20] by the Schwinger mechanism [21] ~ exp(— 7rm^/2K), 
is considerably enhanced. This corresponds either to a 
dramatic enhancement of the string tension k (from the 
default ~ IGeV/fm to 3GeV/fm) or to quark masses 
ruq that are reduced from their constituent quark val- 
ues to current quark values as motivated by chiral sym- 
metry restoration [18]. In a fully analog way, p abun- 
dances can be explained since an enhanced string tension 
similarly leads to an increased production probability of 
antidiquark-diquark pairs which is needed to account for 
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the experimentally observed yields [22]. However, this 
argument, providing additional motivation, has to be re- 
considered if multifusion processes [23], e.g., Stt BB, 
that are neglected in microscopic transport models based 
on 2 — > n scatterings, contribute significantly to the 
baryon pair production. A variation of the string ten- 
sion from K = I GeV/fm to 3 GcV/fm increases the pair 
production probability of strange quarks (compared to 
light quarks) from 7^ = P{ss)/P{qq) = 0.37 to 0.72. 
Similarly, the diquark production probability is enhanced 
from = P{qqqq)/P{qq) = 0.093 to 0.45. In general, 
heavier flavors or diquarks (Q) are suppressed according 
to the Schwinger formula [21] by 

Additional motivation of an effectively enhanced string 

tension in a densely colored environment is given due to 
its relation to the Regge slope a'. Based on a rotational 
string picture the string tension k is related to the Regge 
slope a' by [24,25] 

« = TT-^- (2) 

The empirical value of the Regge slope for baryons is a' ~ 
1 GeV~^ [26] that yields a string tension of approximately 
1 GeV/fm. However, high-energetic processes dominated 
by Pomcron exchange, characterizing the multi-gluon ex- 
change processes existent in high-energetic nucleus nu- 
cleus collisions, are described by a Regge trajectory 
(Pomeron) with a smaller slope of a'p ss 0.4GeV~^ 
[27,28]. According to Eq. (2) this translates into a con- 
siderably larger (effective) string tension k. 

It has been suggested recently that the magnitude of 
a typical field strength at RHIC energies might be as 
large as 5 — 12 GeV/fm [29] (as a result of collective ef- 
fects related to quark-gluon-plasma formation). This in 
turn may enhance an anti-flow component resulting in a 
novel wiggle structure (that in general can have differ- 
ent origins [29-33]) around midrapidity in semi-central 
collisions [29]. 

For our investigation of the baryon dynamics in heavy 
ion collisions at RHIC a microscopic transport approach 
based on the covariant propagation of constituent quarks 
and diquarks accompanied by mesonic and baryonic de- 
grees of freedom is applied [34]. It simulates multiple 
interactions of ingoing and newly produced particles, the 
excitation and fragmentation of color strings and the for- 
mation and decay of hadronic resonances. Standard type 
fragmentation functions are used to assign the individual 
hadrons of the string break-up a longitudinal momentum 
fraction x. Leading nucleons obey the distributions 

f{x)n = exp(-(x - 6)V(2a^)) 

with the maximum at 6 = 0.42 and the width a = 0.275. 
For the produced hadrons a Field-Feynman fragmenta- 
tion function [35] is used 



f{x)p = (1/3)((1 - xfc{d + 1) + 1 - c) 

with constants c = 1 and d = 2. At RHIC ener- 
gies, subliadronic degrees of freedom are of major impor- 
tance. They are treated in the Ultra-relativistic Quan- 
tum Molecular Dynamics (UrQMD) model [34] via the 
introduction of formation times for hadrons that are 
produced in the fragmentation of strings [20]. Leading 
hadrons of the fragmenting strings contain the valence- 
quarks of the original excited hadron. They are allowed 
to interact even during their formation time, with a re- 
duced cross section defined by the additive quark model, 
thus accounting for the original valence quarks contained 
in that hadron [34]. Newly produced (di)quarks are 
not allowed to interact until they have coalesced into 
hadrons. Their formation times are inversely propor- 
tional to the string tension 

tf ~ 1/k. 

The larger the string tension the shorter and short-living 
are the strings with a certain total energy. The cross 
sections for leading quarks and diquarks are set to the 
(later to be formed) hadron hadron cross sections scaled 
down by quark counting 

^g(M)/i = O'M/t/2, <^qq(B)h = '^'^Bh/^, ^'^^'^q{B)h = '^Bh/^, 

where, e.g., q{M) represents a quark of an still unformed 
meson M. In this way the overall stopping behaviour 
of protons in central S-l-S collisions at E\ah = 200 AGeV 
can be described reasonably [36]. The importance of the 
diquark dynamics, i.e., (re) scatterings based on quark 
model cross sections, for the stopping behavior thus be- 
comes apparent [36,37]. These secondary scatterings are 
important for transporting baryon number from pro- 
jectile and target rapidity closer to midrapidity. The 
baryon-baryon collision spectrum is largely dominated by 
diquark degrees of freedom [36]. Only for the very first 
collisions and at the later stage, lower energetic collisions 
involve fully formed baryons. Additional contributions 
arise from (single) quarks carrying baryon number after 
diquark breaking. The diquark breaking mechanism [38] 
is assigned a probability of (~ 10%) . 

Rather exotic mechanisms, e.g., baryon junctions 
[39,40], have been suggested in order to imderstand the 
baryon number transport in high energetic nucleus nu- 
cleus collisions. In this approach, baryon number is not 
necessarily carried by the valence quarks but may reside 
in a non-perturbative configuration of gluon fields. 

There are other important model ingredients to the 
baryon dynamics, i.e., in particular to the baryon trans- 
port over several units in rapidity and the baryon- 
antibaryon pair production. The total and clastic baryon 
antibaryon cross sections at high energies (> 5GeV/c) 
are parameterized according to 

O-tot/eKPlab) =a + + Cln^piab + rflnplab 
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with appropriate constants a, b, c, d, and n to fit exper- 
imental data [34]. The pp annihilation cross section is 
parameterized as 

= Mso/s)[{aho)/{{s - sof + aho) + b]. 

For the pn and the BB cross sections we assume 

and 

with = 1 if the baryons are both nonstrange and 

'^s = o'i^m/'''aqm fo'" baryons with Si > 1 (anti)strange 
quarks, respectively. The energy-independent additive 
quark model cross section in units of mb, (here only used 
to scale the energy-dependent BB cross section), follows 
from 

a^T = 40(2/3)("i+"=)(l-0.4si/(3-ni))(l-0.4s2/(3-n2)) 

with m = 1 for mesons and rij = for baryons. For 
further details of the model we refer to [34] . 

The anti-baryon-to-baryon ratios at midrapidity are 
shown in Figure 1 as a function of the strangeness content 
\S\, both for K = 1 GeV/fm as well as for the strong color 
field scenario (SCF, k = 3GeV/fm). 
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FIG. 1. Antibaryon-to-baryon ratios at midrapidity as a 
function of tlie strangeness content 15*1 in Au+Au collisions 
at RHIC (v^jvjv ~ 200 GeV) . Calculations with a string ten- 
sion of /t = 1 GeV/fm are shown on the left and the results 
with strong color fields (k = 3 GeV/fm) are shown on the 
right. The open and full circles correspond to rather central 
(6 < 3fm) and peripheral (9 < 6 < 12 fm) collisions, respec- 
tively. 

The following systcmatics arise from these calculations, 
(i) The B/B-ratios increase with tlic strangeness content 
\S\ of the baryons. (ii) The i?/i3-ratios increase with 
impact parameter b. There is stronger absorption of an- 
tibaryons in central collisions (reducing the numerator of 



the ratio) and more stopping (increasing the denomina- 
tor), (iii) Most important, the B/B-ratios also increase 
with the color field strength n. The dependence of the 
B / B-T&tios on k is, a priori, not directly transparent. On 
the one hand, an increased pair production drives the ra- 
tio towards 1. On the other hand, reduced formation 
times and, hence, stronger baryon transport, reduce the 
ratios. 

In general, the B/B-r&tio is determined through the 
pair production of BB pairs, i.e., baryons and an- 
tibaryons produced in equal amounts, and the colli- 
sions which transport (net)baryons from beam rapidity to 
midrapidity. Both, the pair production and the number 
of collisions (through smaller formation times) increase 
with increasing field strength k. So in a simple ansatz 
one may write 

B/B^P{x)/[P{x)+C{x)l 

where P{x) is the number of produced BB pairs and 
C{x) is the number of (net)baryons (which were trans- 
ported to midrapidity). Both quantities depend on some 
coupling X, including the color field strength k: (other 
couplings are for example the (re) absorption channels). 
If collisions were neglected, i.e., C{x) = 0, the B/B-ratio 
equals unity. This case corresponds to (freely) decay- 
ing strings and yields the oversimplified transparency as- 
sumption as mentioned in the introduction. However, if 
collisions are taken into account (and this is the case in 
this model where the leading quarks and diquark cross 
sections are always nonzero, even in the formation time,) 
then one obtains B/B-ratios smaller than unity even for 
small field strengths k for which the qq^qq pair produc- 
tion probability is strongly suppressed (see eq. (1)). For 
large field strengths k, pair production dominates result- 
ing in larger ratios B/B ~ P{x)/[P{x) + C(x)] < 1. 

In Table I, the results for the central collisions are com- 
pared to RHIC data from STAR at ^/s = 130 GeV (10% 
most central collisions) [41]. The results from other col- 
laborations agree within error bars. Obviously, the calcu- 
lations obtained with strong color fields agree better with 
the data. Note that feeddown corrections, e.g., protons 
from A decays, may further reduce the experimentally 
(uncorrected) ratios. Experimental error bars are statis- 
tically only. 

TABLE 1. p/p, A/A, E/E, and fi/fi ratios at midrapidity 
for central Au-|-Au collisions at RHIC for calculations with 

n — 1 GcV/fni. with strong color fields and compared to (par- 



tially preliminary) STAR data at ^/sj^ 



130 GeV. 





K = 1 GeV/fm 


(SCF) /t = 3GeV/fm 


exp. data 


p/p 


0.41 ± 0.01 


0.52 ± 0.01 


0.72 ± 0.01 


A/ A 


0.56 ±0.01 


0.66 ± 0.02 


0.73 ± 0.03 




0.64 ±0.03 


0.72 ±0.04 


0.83 ±0.03 


n/fi 


0.60 ±0.12 


0.80 ± 0.08 


0.95 ±0.15 
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The rapidity spectra for netprotons {p — p) and net- 
baryons {B — B) are shown for different centrahties in 
Figure 2. For the most central collisions, the rather 
strong maxima at y ~ 3.2 (k = IGeV/fm) are weak- 
ened and shifted closer to midrapidity {y ~ 1.8) in the 
case of strong color fields (k = 3GeV/fm). The num- 
ber of netprotons at midrapidity seems to be indepen- 
dent on the color field strength 12) while the number 
of netbaryons increases (from about 32 to 42). Hence, 
the larger netbaryon number at midrapidity results from 
strange baryons. 
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ues [42] too small compared to data at RHIC [43]. 

The transport of other quantum numbers that can 
be examined in AA as well as in pA collisions include 
strangeness [18] and isospin [34,44]. Those can provide a 
complementary picture of the baryon number transport 
and do not need to be conserved locally. 

We have calculated the antibaryon-to-baryon ra- 
tios, netbaryon and netproton rapidity distributions for 
Au+Au collisions at RHIC energies, and have demon- 
strated the important effects of strong color fields on 
these observables. An enhanced pair production of 
strange quarks and diquarks and reduced formation times 
lead to larger B/B-iatios and to a stronger rapidity loss 
for the netprotons and netbaryons. More netbaryons 
are transported to midrapidity in the strong color field 
scenario which is accompanied by a hyperonization of 
these netbaryons. This hyperonization can be observed 
through an enhanced (A — A)/(p — p) ratio and possibly 
represents the clearest signature of the flavor equilibra- 
tion in relativistic heavy ion collisions. 
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FIG. 2. Rapidity distributions for nct-baryons (top) 
and net-protons (bottom) in Au+Au collisions at RHIC 
(\/*jvAr = 200 GeV). Calculations with a standard string 
tension (k = IGeV/fm) are shown on the left and the re- 
sults with strong color fields {k = 3GeV/fm) are shown on 
the right. The symbols correspond to different centralities 
(6 < 3fm (circles), 3 < fe < 6fm (squares), 6 < 6 < 9fm 
(triangles), and 9 < fe < 12 fm (diamonds)). 

This hyperonization can also be seen in the (A— A)/(p— 
p) ratio at midrapidity which is ^ 0.45 for calculations 
with K = 1 GeV/fm and strongly enhanced ~ 0.8 for the 
strong color field (SCF) (k = 3GcV/fm) scenario, corre- 
sponding to the experimental value [4,1]. Note that S^'s 
and T,^'s are included in the A's and A's, respectively, 
because they cannot be distinguished experimentally. 

Another consequence of strong color fields is an in- 
crease of the pressure (gradients). This, in turn, may 
lead to strongly enhanced elliptic flow values as large as 
observed in the data. This provides a possible explana- 
tion for the lack between calculations that assume the 
vacuum string tension and that predict elliptic flow val- 
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